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ABSTRACT 

The mechanism for gamma ray bursters and the detection of gravitational waves 
(GWs) are two outstanding problems facing modern physics. Many models of gamma 
ray bursters predict copious GW emission, so the assumption of an association between 
GWs and GRBs may be testable with existing bar GW detector data. We consider 
Weber bar data streams in the vicinity of known GRB times and present calculations 
of the expected signal after co-addition of 1000 GW/GRBs that have been shifted to 
a common zero time. Our calculations are based on assumptions concerning the GW 
spectrum and the redshift distribution of GW/GRB sources which are consistent with 
current GW/GRB models. We discuss further possibilities of GW detection associated 
with GRBs in light of future bar detector improvements and suggest that co-addition 
of data from several improved bar detectors may result in detection of GWs (if the 
GW/GRB assumption is correct) on a time scale comparable with the LIGO projects. 
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1 INTRODUCTION 



The attempts to detect gravitational radiation began with 
the work of Weber (1960). The existence of gravitational 
waves has been verified indirectly by observations of the de- 
cay in orbital period of the pulsar PSR 1913-1-16 (Hulse & 
Taylor 1975) , for which the 1993 Nobel Prize in Physics was 
awarded. However, direct detection continues to prove elu- 
sive even though some authors believed it to be possible 
before the turn of the century (eg. Thorne 1992). Detection 
remains difficult due to the nature of the coupling between 
energy and spacetime. Although a large amount of energy 
might be stored in gravitational waves (GWs), the corre- 
sponding amplitude of the waves, measured in dimensionless 
strain, h, is exceedingly small due to the very small magni- 
tude of the proportionality constant in Einstein's equations: 
a large stress-energy produces a very small metric perturba- 
tion. 

The most sensitive detectors currently in operation are 
of the Weber bar variety. The main detectors are maintained 
at the University of Western Australia (eg. Heng et al. 1996), 
Louisiana State University (eg. Mauceli et al. 1996), the 
University of Rome (eg. Astone et al. 1993), Legnaro Na- 
tional Laboratories (eg. Cerdonio et al. 1995) and Labora- 
tori Nazionali di Frascati (eg. Coccia et al. 1995). These 
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detectors have a very low bandwidth as they rely on detec- 
tion of the resonance in the bar (or sphere) after a burst of 
GWs. Near the resonance frequency, bar detectors are ca- 
pable of reaching sensitivities of ~ 10~^^-10~^^ (eg. Tobar 
et al. 1995) and future sensitivities are expected to reach 
~ 10"^° (eg. Tobar et al. 1998). Spurious, non-thermal ex- 
citations of resonant bars typically number ~ 40 per day 
(Heng et al. 1996). Thus, detection of a real, individual 
GW burst will require coincidence with detections by other 
instruments around the world. No convincing coincidences 
have been recorded to date. 

The LIGO/ VIRGO network is expected to detect grav- 
itational waves (Finn & Chernoff 1993; Cutler & Flanagan 
1994; Thorne 1994; Flanagan & Hughes 1998), if not in the 
initial stages when the sensitivity of LIGO is ~ 10^^^, then 
in the advanced stage when it is ~ 3 x 10^^''. This will mark 
the beginning of GW astronomy, promising views of neutron 
star (NS) and black hole (BH) coalescence and collisions. It 
will also yield detailed information about the dynamics of 
pulsar systems and may even allow us to study infiation and 
the Big Bang itself more directly and in much more detail 
than other astrophysical techniques allow. 

Another problem that has confronted science for 
decades has been the origin of gamma ray bursts (GRBs). 
Many theories for the emission of GRBs centre on a mech- 
anism involving large mass quadrupole moments. Hence a 
wide variety of GRB models predict large GW luminosities. 
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For example, NS mergers are good candidates for GRBs (eg. 
Kochanek & Piran 1993; Dokuchaev, Eroshenko & Ozernoy 
1998) and are expected to result in ~ 10^^ ergs in GWs 
(Centrella & McMillan 1993; Ruffert & Janka 1998). Such 
energies result in detectable GW signals if the NSs are in 
the Galaxy but if the source is cosmological both modern 
bar detectors and the advanced LIGO system cannot detect 
such events. 

It is this expected association between the two phenom- 
ena that has prompted the present investigation. With the 
expectation that GW bursts are temporally correlated with 
GRBs, we may co-add existing GW data surrounding mea- 
sured GRB events. This would result in a composite GW 
data stream with all GRB events shifted to a common 'GRB 
time'. We explore this idea below and present a method of- 
fering up to a factor of ~ 30 increase in signal-to-noise ra- 
tio (SNR) with the co-addition of 1000 existing GW data 
streams surrounding measured GRB times. Throughout the 
work, we assume that all the measured GRBs are associated 
with GW emission. 

In this article we calculate the expected signal from a 
typical bar detector assuming that the bar is noiseless. Our 
aim is to show directly, with some generality, whether avail- 
able data from GRB and GW detectors can be used to detect 
GWs before the first LIGO/ VIRGO detection. Such a detec- 
tion would also be confirmation of an association between 
the GRB mechanism and GW emission. 

Our calculations are based on various assumptions con- 
cerning source distribution in spacetime, the nature of the 
source and the corresponding GW emission, all in line with 
current models. We compare our calculations with the noise 
floor of the bar and decide on the viability of such a method 
for detection using bar detectors. We also discuss the viabil- 
ity of detection by the LIGO system. 

The paper is organized as follows. In Section 2 we pro- 
vide the formalism for our calculations, discussing the as- 
sumptions concerning the GW source and detector. Vari- 
ous models are considered in Section 3 that are consistent 
with the assumption that GRBs are associated with GW 
emission. From this discussion we extract features such as 
the redshift distribution of sources and the spectrum of the 
source. In Section 4 we compute the signal seen at the detec- 
tor for a single GW/GRB. We also simulate the co-addition 
of bar data with the assumption of a common temporal zero 
point and we discuss the results for the various models in 
Section 4. In Section 5 we discuss the relevance of our result 
for the LIGO/ VIRGO experiments and provide our conclu- 



detector we need to find the GW spectrum at the bar due 
to a source under this distribution. Therefore, let us assume 
that the frequency of GWs at the source is described by the 
spectral profile f[u)dv. Again, J f[u)dv — 1. Also, since we 
are considering the spectrum of the source, we may calibrate 
the GW flux scale by assuming that the average total en- 
ergy emitted in GWs by the source is F. We explain our 
reasoning below. 



2.1 GW emission and propagation 

Our first assumption is of a distribution of mean total GW 
fluxes at the detector. The analysis is based on the as- 
sumption that GRBs are associated with GW emission. We 
further assume that the observed luminosity distribution of 
GWs in the detector frame follows that of GRBs. This as- 
sumption is valid only if the mechanism for gamma ray and 
GW emission is closely physically related. We make this as- 
sumption on the basis that our results may then be used 
as constraints for such a subset of GW/GRB theories and 
because of the lack of other reasonable models for the lumi- 
nosity distribution of GW sources under the GW/GRB as- 
sumption. The knowledge of whether this relationship holds 
in reality would be invaluable for theories of gamma ray 
bursters. 

It is common to define the Earth frame luminosity dis- 
tribution of GRBs in terms of the number of bursts with 
peak gamma ray flux greater than a given peak flux, P, 
as Nc,rb{> P)- The fact that this distribution strays from 
a —3/2 power law at low P is an indication of the cosmo- 
logical origin of GRBs. This has been conflrmed by several 
optical transient observations and redshift determinations 
(eg. Metzger et al. 1997; van Paradijs et al. 1997; Djorgov- 
ski et al. 1998; Kulkarni et al. 1998). P is usually defined in 
terms of a given Earth frame energy band and the integra- 
tion time over which the flux was allowed to accumulate. We 
use the distribution measured by the BATSE instrument - 
supplied in the BATSE 4Br catalogue (Paciesas et al. 1999) 
- choosing the 50-300 keV energy band and the 64 ms inte- 
gration time. In our calculations, it is convenient to flnd the 
cumulative distribution of total GW fluxes. We deflne this 



n™.[< P(5')] = l-^^^=^#^ 

JVgrb 



(1) 



where A^grb is the total number of GRBs in the distribution 
and P{S') is defined as 



2 FORMALISM 

Our aim is to calculate the signal seen at the output of a 
GW bar detector, with a single resonance band, due to a 
single GW burst. However, we keep in mind that we ulti- 
mately wish to simulate the composite effect due to many 
bursts (shifted to a common burst time) ; we wish to obtain 
a Monte Carlo estimate of the signal at the bar output af- 
ter data processing. We therefore assume that there exists 
a distribution of GW/GRB sources in spacetime according 
to the normalized redshift distribution, (j){z)dz. That is, we 
require that J (f>{z)dz — 1. To calculate the effect on the 



Pis') ^j;S'. 



(2) 



Here, P is the average value of the peak flux in gamma rays 
and S' is the mean flux in GWs at the detector. 

P is easily computed from the data in the BATSE 4Br 
catalogue. We obtain 5*' by evolving F with cosmology from 
the source to the detector. This is particularly simple in the 
case of a flat universe in which the cosmological constant is 
vanishingly small (Misner, Thorne & Wheeler 1973, p. 783), 



F 



where the distance, ^^.(zni) is deflned by 



(3) 
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7^(Zm) = 



2 + 2Zm - 2V1 + Zn 



(4) 



Here, //o is the Hubble constant and Zm is the mean of the 
redshift distribution (c has been set to unity in Eqs. |^ and 
Throughout this work, we use Hq = 75 kms~^Mpc~^. 
So, by choosing a model of the source, we can estimate F 
and normalize the flux scale of ricw (< P{S')) to the units of 
GW flux using Eq. |[ 

Since we are interested in the frequency space at the 
detector, we may convert the assumed redshift distribution, 
(j){z)dz, to a GW frequency distribution at the detector, 
(l)'^{v')dv' , using the simple transformation, 1 + z = v jv' . 
This leads directly to 



(5) 



provided that we have a value for v, the frequency of GW at 
the source at redshift z. We interpret v as being the mean 
frequency at the source under the distribution j(y)dv. It is 
also convenient to define the cumulative frequency distribu- 
tion at the detector, 



(6) 



Since we ultimately wish to co-add the signals from 
many GWs bursts associated with GRBs, we may choose 
a luminosity from nGw(< S') by randomly selecting a value 
for tigw and finding the corresponding value of S' . Similarly, 
we use Eq. ^ to find a random value for the mean GW fre- 
quency at the detector. Since we have assumed a normalized 
GW spectral profile of the source, f(u)dv, which has a mean 
frequency of u, we may find the mean GW frequency at the 
detector given a single GW/GRB burst. That is, by speci- 
fying both the mean frequency at the source and at the de- 
tector, we have randomly selected the source redshift. Using 
this redshift we can find the full spectrum at the detector 
by multiplying our selected value for S' by the redshifted 
spectral profile, f^{i'')di/' . 



2.2 GW detection 

2.2.1 General 

We may now discuss the detection of such a GW burst. We 
assume that the GW burst bandwidth is much greater than 
that of the bar detector since the later is typically of the 
order of ~ 0.3 Hz (Tobar et al. 1995). The cross section for 
a single bar resonance in response to GWs in such a situation 
is given by (Paik & Wagoner 1976) 



32M /G 
15 



_ f—\ (- 

TV \ C J \C 



l + O 



(7) 



where M is the mass of the bar, v is the speed of sound in 
the bar, R is the radius of the physical cross section and L 
is the length of the bar. Eq. ^ assumes that, on the average, 
GW bursts are incident isotropically on the detector with 
random polarisations. In most bar detectors, R <^ L and 
so the second order terms are neglected. We can then find 
the amplitude of vibration (in strain) that the absorption of 
GWs causes in the bar (Misner, Thorne & Wheeler 1973, p. 
1039), 



-I 



S'^,,Jl{u')dv' 
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(8) 



2.2.2 The Niobium bar detector at UWA 

The discussion above has remained general for most cylindri- 
cal bar detectors. However, for the remainder of the analysis, 
we take the specific example of the Niobium bar detector at 
the University of Western Australia, NIOBE, since we wish 
to model the detection algorithms in order to estimate the 
signal seen at the outputs. The results will carry a factor of 
~ 5-10 uncertainty and so they will remain representative 
of most bar detectors. In particular, we consider the sim- 
plest algorithm to analyse GW data streams, the zero order 
prediction (ZOP) technique. Here we summarize some main 
points. For details of the detection methods used for NIOBE 
and for descriptions of the ZOP algorithm, see Tobar et al. 

(1995) , Dhurandhar, Blair & Costa (1996) and Heng et al. 

(1996) . 

GW data are sampled at 0.1 s intervals and passed 
through an optimal low-pass filter before being decimated 
to 1 Hz samples. Correction factors for the low-pass proce- 
dure are included and the results are scaled to noise tem- 
perature or strain. For an ~ 1 ms burst of GWs, the rise 
time for the amplitude of vibration in the detector is also 
~ 1 ms. However, due to the high Q-factor, the amplitude 
decays exponentially over a time scale of ~ 200 s. The fact 
that the data are decimated means that the signal seen after 
the data are processed is essentially the time derivative of 
the data. That is, if we model the decay in strain amplitude 
by 

I / dt \ f dt 

exp I —t— — exp I —t- 



h{t) = V2h 



Tin 



(9) 



then we may obtain an estimate of the signal after data 
processing using 

dt ( dt 



h(t) = V2h 



T 



■ exp 
dt 



\ mode / 

dt 



dt ( dt 
exp I —t 

Tint V Tint 



(10) 



where f is a dimensionless measure of the time after the 
GW burst that takes into account the sampling time in the 
decimated data. Tmode is the ring down time for the mode 
of oscillation in the bar. Tint is the integration time for data 
samples and dt is the time step for data sampling. In our 
case dt — Tint and Tmodo ~ 200 s. This analysis does not 
take account of the noise floor of the detector which will 
completely swamp such a single GW burst signal. Also, we 
do not take account of any attenuation of the signal resulting 
from the low-pass filtering; we assume that the processing 
algorithm corrects for this. 

2.3 Co-addition Simulation 

We are now in a position to set up a Monte Carlo simulation 
to estimate the signal in Eq. ^| We must choose a model 
for the source so as to estimate the GW spectral profile and 
the mean total GW energy output per burst. The formalism 
above implicitly assumes that all GW bursts last for less 
than ~ 1 s and this must be reflected in the model of the 
GW emission. We must also assume a redshift distribution 
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of sources. The models considered here wiU be discussed in 
Section 3. Given one of the set of models, however, we may 
find a value for S' fl{v')dv'\^i and use Eqs. ]A and H to find 
the signal after the data processing, Eq. hcT By using the 
cumulative distributions, Eqs. |l] and ^, we can build up a 
simulation of many bursts and find their collective effect on 
the detector. 

Throughout this work we assume, as a worst case esti- 
mate, that GRB times are uniformly distributed over 10 s af- 
ter the emission of the corresponding GW bursts (Kochanek 
& Piran 1993). However, we provide a scaling of our results 
for a range of time lags in Section 4 (see Fig. 2) since a less 
conservative view would put most GW bursts only ~ 1 s 
before the GRB. The 'GRB time' may be the trigger time 
in the BATSE catalogue or it may be the time at which 
the gamma ray flux sees a signiflcant increase over thresh- 
old. The difference between such deflnitions could be, at 
most, ~ 10 s. This issue is complicated by the irregular and 
complex intensity proflles of GRBs (eg. Fishman & Meegan 
1995). Here, for purposes of simulation, we assume that this 
problem is resolved so that all the 'GRB times' are aligned. 

Finally, we note that the treatment above is valid only 
for a bar detector with a single resonance band. However, 
as our calculations are aimed at simulating the effect on 
NIOBE, which has two resonant frequencies, the plus mode 
at 1^0+ ~ 713 Hz and the minus mode at Vq_ ~ 694.5 Hz, 
we make the following adjustment to the formalism. NIOBE 
acts approximately as a coupled two mass oscillator. As a 
rough approximation, we may therefore model the detector 
as having two resonance bands with cross sections given by 
Eq. H centred at the relevant resonant frequencies. 



3 MODELS 

In this section we consider the parameters associated with 
various GW/GRB models. Most of the resulting consistent 
sets of parameters describe binary coalescence models as 
these seem to be the most likely GW/GRB source if indeed 
the GW/GRB assumption is correct (Kochanek & Piran 
1993). The remaining sets are, therefore, more speculative 
and provide useful comparisons with the binary coalescence 
models in view of the flnal results in Section 4. 

In order to define a specific GW/GRB model we must 
specify the following. 

(i) Burst Mechanism: Among the competing mecha- 
nisms for GRBs are NS collision and coalescence (e.g. Piran 
1991; Kochanek & Piran 1993; Totani 1997; Ruffert & Janka 
1998) and other more general coalescence scenarios, such as 
BH-NS (eg. Lee & Kluzniak 1998) and white dwarf-BH co- 
alescences (eg. Fryer et al. 1999). These mechanisms all im- 
ply that the source space density should follow (if somewhat 
lagged in time) the massive star formation rate (eg. Totani 
1997; Krumholz, Thorsett & Harrison 1998) . Other authors 
prefer to assume source space densities that are similar to 
the quasar redshift distribution (eg. Horak, Emslie & Hart- 
mann 1995) without venturing as to what the mechanism 
for GRBs might be. 

(ii) Spectral profile, f{v)dv: For close binary systems, 
the emission of GWs leads to energy loss and in-spiral. The 
GW frequency therefore follows a 'chirrup' pattern, mono- 
tonically increasing with time. The Newtonian approxima- 



tion for the in-spiral is a good one up until the few rotation 
periods before coalescence when post-Newtonian corrections 
become important. The innermost stable orbit marks the 
point at which GW emission begins to decrease because of 
the increase in the infall rate of the component stars. The 
innermost stable orbit occurs at an orbital separation of 
r ~ 6M (Kidder, WiU & Wiseman 1992, setting G = c = 1) 
and so Cutler & Flanagan (1994) propose that the GW emis- 
sion will shut off at roughly u = (6'^''^7rA/)^^ for M the to- 
tal mass of the binary system. The important parameters 
to obtain here are the mean frequency of GW emission at 
the source and the spectral profile. For a NS-NS system, the 
mass may lie in the range M = (2.70 ± 0.08)Mq (Thorsett 
& Chakrabarty 1999). We therefore find that the frequency 
of maximum emission is i/ ~ 1600 Hz and that the width 
of the distribution of these frequencies is ~ 100 Hz. Assum- 
ing that the 'chirrup' is very steep near the coalescence and 
that the frequency of emission is within ~ 500 Hz of v for 
^ 1 s (eg. Allen 1997; Allen 1999) , we may assume, as we 
ultimately wish to co-add signals at the detector, that the 
spectral profile at the source for NS-NS coalescence has a 
width of cr ~ 500 Hz. 

The issue of the actual value of v is of more importance 
here since we are dealing with such a low bandwidth detec- 
tor. However, it is a poorly known quantity due to the slow 
convergence of post Newtonian corrections to the estimate 
above. Using the results of post Newtonian calculations of 
the final moments of in-spiral (Finn & Chernoff 1993; Cut- 
ler & Flanagan 1994; Blanchet 1996) we find values in the 
range ~ 1200 to ~ 1500 Hz. We also note that estimates 
by Kidder, WiU & Wiseman (1992, 1993), Damour, Iyer & 
Sathyaprakash (1998) and Allen et al. (1999) suggest that 
V = 1420, 2042.6 and 1822 Hz respectively for M = 2.8Mq. 
It is therefore desirable to assume a range of frequencies for 
each model of binary coalescence. 

We also note that, for coalescence scenarios involving a 
BH, the higher mass of the BH relative to that of a NS 
reduces the GW frequency. It is therefore likely that GW 
emission from such systems will result in a low flux at the 
detector resonant frequency. Consequently, we do not con- 
sider such systems in our calculations. 

With regards to NS collisions, little information can be 
gained as to the general spectral profile of the source since 
this depends on the separation of the path trajectories. 
Therefore, we do not consider collisions in our calculations. 

Finally, assuming that the GW/GRBs are associated with 
the quasar redshift distribution gives no information on the 
distribution of frequencies at the source; the GW mechanism 
has not been specified. We will therefore assume a broad 
band emission in such an instance so as to compare with re- 
sults from the NS coalescence models. We will also consider 
a range of mean redshifts, Zm, in such models. 

(iii) Mean total emitted GW energy, F: For NS coa- 
lescences, this is known to reasonable accuracy, but clearly, 
for our quasar model, we have no information on this pa- 
rameter. Let us therefore remove F from consideration in 
this instance by assuming that it is the same as the range 
we will assume for NS coalescences. Various authors sug- 
gest that F lies in the range 10^^-10^"^ ergs (eg. Rasio & 
Shapiro 1992; Centrella & McMillan 1993; Oohara & Naka- 
mura 1995; Lee & Kluzniak 1998) based on both numerical 
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and approximate theoretical grounds. We therefore provide 
calculations for this order of magnitude. 

(iv) Number of bursts, N: From the formalism in Sec- 
tion 2, we see that if the number of bursts is large then the 
sum of all h for all bursts will scale approximately linearly 
with N. Therefore, throughout our calculations we will as- 
sume that A'^ = 1000 as a representative value. We compare 
this with the number of bursts in the BATSE 4Br cata- 
logue, A^BATSE = 1637, and note that over 2000 burst times 
have been recorded. If the co-addition of GW data were to 
be performed, the number of bursts available would be re- 
stricted by the fact that most bursts do not have a simple 
intensity profile over the BATSE integration time. That is, 
it may be difficult to assign a burst time to a given GRB 
due to the complex nature of the profile. We therefore keep 
A'' = 1000 as an estimate of the number of useful bursts for 
the current purpose. 

In any co-addition of "real" data, the SNR will rise in 
proportion to \/iV for A'^ the number of data streams to be 
co-added. The noise is co-added together with the signal; 
the signal rising in proportion to A'^ and the noise in pro- 
portion to \/iV. However, our co-addition of simulated data 
does not include the noise characteristics of the bar. If we 
are to compare our final results with the sensitivity of the 
detector for a single burst, the comparison must be between 
our noiseless results and y/N times the typical noise at the 
resonance. We make this comparison in Section 4. 

We have therefore narrowed our view of models of 
GW/GRBs to two main sets of parameters. The first is the 
set associated with the assumption that GRBs are due to 
coalescence of binary NS systems. We shall call this the 
NS-NS model and take those parameters discussed above 
for it. We assume a star formation redshift distribution that 
follows that of Madau, Pozzetti & Dickinson (1998) (here- 
after MPD) and compare this with that given by Pascarelle, 
Lanzetta & Fernandez-Soto (1998) (hereafter PLFS). The 
later compares well with the findings of Totani (1997) that 
the star formation rate needs to be a factor of 5-10 higher 
at redshifts 2 > 2 for the NS-NS scenario to be consistent 
with the BATSE data. 

However, as a useful comparison, we parametrize the 
somewhat undefined quasar model by assuming a Gaussian 
redshift distribution, centred on a range of redshifts, 1.5 < 
< 3.0 with a constant width, a — 0.45 (Horack, Emslie 
and Hartmann 1995). The range in Zn-i may be interpreted 
as a time lag factor (positive or negative) due to the specific 
(unknown) mechanism of GW/GRBs. We also include such 
a time lag eff'ect in the star formation redshift distributions 
due to the length of time needed for complete in-spiral of 
close NS-NS systems. We vary this time lag in redshift space 
from ~ 0.8 to by allowing the redshift at which the star 
formation rate first reaches a maximum, Zp, to vary. We 
illustrate some typical redshift distributions in Fig. 1. 



4 RESULTS AND DISCUSSION 

For a specific model defined in the previous section, with a 
given mean source frame frequency, v, and a specific redshift 
distribution, we obtain results typified by those shown in 
Fig. 2. Of note is the distribution of h shown by the points. 
These are co-added to give a generally increasing signal. 




Redshift z 

Figure 1. Comparison of three typical redshift distributions. The 
top panel shows a MPD-like distribution and the middle panel 
shows a rough fit to the redshift distribution in PLFS. The bottom 
panel shows the quasar-like distribution used as a comparative 
example in our calculations. All distributions are normalized to 
unity. 

h{t), up until the common GRB time. At times t > 0, the 
collective derivative patterns (Eq. |l^) of all the events cause 
the decay of the signal. The inset of Fig. 2 shows the scaling 
for different values of the time lag between the GW burst and 
the GRB. Although we have calculated this scaling with the 
same model parameters as those described in the caption to 
Fig. 2, the scaling is model independent and can be applied 
to any of our subsequent results (Figs. 3-6). 

Due to the fast decay time (~ 1 s) of the signal out- 
put from the ZOP algorithm (Eq. |l^, the time lag param- 
eter is clearly an important one. However, in a real data 
co-addition, the separation between the GW burst time and 
the common GRB time is also subject to other uncertain- 
ties as described in Section 2. We therefore prefer to assume 
that these uncertainties limit us to an effective maximum 
difference between the GW and GRB times of 10 s and we 
use this value in calculating all subsequent results. 

So as to give information for all models as functions of 
and the variable redshift parameter (zp in the case of the 
NS-NS models and Zm in the case of the quasar model), our 
results are presented as three-surfaces in parameter space 
in Figs. 3-6. Fig. 3 shows the results for the NS-NS model 
with a MPD redshift distribution. The vertical axis is the 
peak measured strain in the detector, max(/i(t)), in units of 
10"^"^. We note the sharp rise in signal with increasing mean 
source frequency and the broad leveling off for frequencies 
as high as ~ 1800-2200 Hz. We also see a general increase 
in peak signal with decreasing peak redshift which becomes 
less important at high source frequency. 

The noise on the surface is due to the fact that we 
only use 1000 bursts in the calculations. Of course, we could 
have increased the number of bursts to reduce the fractional 
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Figure 2. Demonstration of the co-addition of 1000 GW bursts. 
We used the NS-NS scenario with the star formation redshift 
distribution of MPD (top panel of Fig. 1) with u = 1600 Hz and 
F = 10®^ ergs. The plot shows each burst as a point on the h scale 
and the integrated signal seen after data processing is shown as 
the solid line. We have assumed that burst times are distributed 
uniformly in time for the 10 s before the common GRB time. The 
inset shows the maximum values of h{t) for a range of time lag 
values relative to that for the 10 s lag, r(lOs). For clarity, we 
have averaged the maximum strain over 10 simulations to obtain 

r(ios). 

fluctuation after the surfaces had been normalized to 1000 
bursts. However, the noise is illustrative of the expected fluc- 
tuations in signal due to the various distributions (spacetime 
distribution, distribution in frequency space etc.) that we 
assume the sources follow. We do note, however, that the 
expected error in our calculations may be as high as a fac- 
tor of 5-10 due to our basic treatment in Sections 2 and 3 
(Misner, Thorne & Wheeler 1973, pp. 1004-1044). That is, 
we have not considered such parameters as the bandwidth of 
the detector, aberrations due to coupling of the displacement 
sensors to the bar etc. For this reason, we shall compare our 
results with the total noise in the NIOBE detector later in 
this section. 

In Fig. 4 we show results for a similar set of parameters 
(see Section 3) but we use the PLFS redshift distribution. 
We note the immediate qualitative difference in that the re- 
sults show less variation in Zp space and that there is a linear 
increase in the signal with increasing frequency. This is sim- 
ilar to the almost linear increase we see in Fig. 3 (for the 
MPD model) at low frequencies and is due to the fact that 
more sources lie at higher redshift in the PLFS model. That 
is, the frequency of the leveling off will be higher since the 
GWs, on the average, at an average redshift in the PLFS 
model, are redshifted below the resonance band in the de- 
tector. Quantitatively though, the PLFS model provides a 
lower average Hux in the resonance band, and we see that 




Figure 3. Three-surface showing the peak measured signal at the 
detector outputs (after data processing), max(/i(t)), as a function 
of the peak redshift, Zp, and the mean source frame frequency, P, 
for the NS-NS model with MPD redshift distribution. Model pa- 
rameters are described in Section 3. We used a mean total emitted 
GW energy at the source of F = 10^^ ergs but the vertical axis 
can be scaled as \/~F for generalization to other values. 

the signal is ~ 1/2 that given by the MPD model throughout 
the parameter space. 

Finally, in Figs. 5 and 6, we provide the results for the 
quasar model. For Fig. 5 we used a uniform frequency dis- 
tribution at the source with an order of magnitude width. 
All other parameters are kept the same as those used for 
the NS-NS models with the exception of the range of red- 
shift parameter used. Fig. 6 shows a similar calculation us- 
ing a Gaussian frequency distribution with the same width 
(a = 500 Hz) as those used in the NS-NS models. We see 
a much slower variation of the signal over the parameter 
space, particularly at lower frequencies. Quantitatively, we 
see that the quasar models have comparative signal level 
with the star formation models and so we find that the sig- 
nal seen at the bar detector is reasonably insensitive to the 
redshift and source model. 

We may compare the results in Figs. 2-6 with the sensi- 
tivity curves provided by Tobar et al. (1995) for the NIOBE 
bar detector. The sensitivity curves show that, within the 
two resonance bands, the root-mean-square noise reaches a 
level of ~ 2 and ~ 7 x lO"'^", for the plus and minus mode 
respectively. For 1000 co-added bursts, the noise rises by a 
factor of VlOOO; ~ 6 and ~ 22 x 10" ^'^ in each band respec- 
tively. However, the maximum signals calculated in Figs. 3-6 
only range from 16-31 x 10"^''. That is, the maximum ex- 
pected signal-to-noise ratio after co-addition is ~ 3 x 10~*. 
Therefore, despite the potential gain in signal evident in the 
co-addition results above, the NIOBE detector and other 
bar detectors, will still be unable to detect the signal with 
1000 GRB times and existing GW data. 

However, plans exist for improvements to be made to 
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Figure 4. Three-surface showing the peak measured signal at the 
detector outputs, max(/i(4)), as a function of the peak redshift, 
Zp, and the mean source frame frequency, P, for the NS— NS model 
with PLFS redshift distribution. Model parameters are described 
in Section 3. We used a mean total emitted GW energy at the 
source of _F = lO^'^ ergs but the vertical axis can be scaled as \/F 
for generalization to other values. 

the NIOBE detector. A 2-mode sapphire transducer is to 
replace the existing vibration sensors and the cryogenic am- 
phfier is to be significantly improved. In addition, the (mi- 
crowave source) pump oscillator driving the transducer is to 
be improved. Tobar et al. (1998) report that these changes 
will decrease the detector noise to ~ 5 x 10"'^^ in strain. A 
more important parameter will be the increase in the band- 
width of the detector from ~ 0.3 Hz in two bands to a single 
band with a width of ~ 50 Hz. We discuss the prospects for 
the improved detector in Section 5. 



5 CONCLUSIONS 

Using models of the GW/GRB source and the spacetime 
distribution of sources, we have estimated the signal seen in 
the processed data of a typical bar detector, NIOBE. As- 
suming a knowledge of the time of the GW bursts (in this 
case, by assuming an association between GW and GRBs) 
allows us to co-add simulated GW data from many bursts. 
Figs. 3-6 show the variation of the resultant signal with the 
main source parameter, the mean GW frequency v, and the 
main spacetime parameter, the peak or mean redshift, Zp or 
Zm - These results assume co-addition of 1000 measured burst 
events. The largest signal that can be expected occurs within 
the NS-NS model with the MPD redshift distribution. The 
peak in signal occurs over a large region of parameter space 
covering mean rest frame GW frequencies v > 1600 Hz and 
all realistic peak redshifts Zp. The maximum expected SNR 
is ~ 5 X 10~^ in this case. Thus, comparing the signal with 
the noise floor of NIOBE, we find that the co-addition, with 




Figure 5. Three-surface showing the peak measured signal at the 
detector outputs, max(/i(t)), as a function of the mean redshift, 
Zm, and the mean source frame frequency, P, for the quasar model. 
Model parameters are described in the text and are set so that 
we may compare these results with those in Figs. 3 and 4. In 
particular, we used a mean total emitted GW energy at the source 
of _F = 10®^ ergs but the vertical axis can be scaled as y/F for 
generalization to other values. However, here we have assumed a 
uniform frequency profile at the source with a width of an order 
of magnitude. 

the present sensitivity of NIOBE and the present number of 
bursts, is insufficient to allow detection by a method based 
on the formalism in Section 2. 

However, the improvements to NIOBE will result in an 
order of magnitude sensitivity gain. The bandwidth will also 
increase by a factor 50.0/(2x0.3) = 83. Tobar & Blair (1995) 
and Tobar (1997) report that, to a good approximation, h oc 
VAf for A/ the bandwidth of the detector. That is, from 
Eq. |8|, a linear change in the cross section of the bar, E^/^, 
corresponds to a change in /i oc y'E^/ . We can therefore find 
an optimistic estimate of the number of bursts, A^b, required 
for us to achieve a SNR ~ 1, 

1000 (^)' (3^)' (^)' 

^10ff_crgs^ ^ 0.6 Hz -^ 

The first term reflects the fact that we have achieved a SNR 
~ 5 X 10~* using 1000 co-added bursts. The second term 
takes account of an increase in sensitivity; Snew is the strain 
sensitivity in the resonance band of the improved detector. 
The flnal terms incorporate changes in the time lag factor 
introduced in Fig. 2, an increased average energy output in 
GWs at the source and a change in the bandwidth of the 
detector respectively. Thus, the number of GW/GRBs to be 
observed could be as low as ~ 8000-10000 in the best case^ 

t That is, we use F ^ 10^^ ergs and a time lag of 2 s between GW 



© 1999 RAS, MNRAS 000, 0-| 



8 M. T. Murphy, J. K. Webb and I. S. Heng 




Figure 6. Three-surface showing the peak measured signal at the 
detector outputs, max(h(t)), as a function of the mean redshift, 
Zm, and the mean source frame frequency, P, for the quasar model. 
Model parameters are described in the text and are set so that we 
may compare these results with those in Figs. 3 and 4. Note that, 
in this case, we also keep the frequency profile as a Gaussian with 
a = 500 Hz at the source. Also, we used a mean total emitted 
GW energy at the source of F = lO^'^ ergs but the vertical axis 
can be scaled as \/~F for generalization to other values. 

Under this assumption, detection of GWs by NIOBE could 
be found in ~ 25 years as the GRB detection rate is ~ 1 
per day. A disadvantage here is that the accumulation of 
new observed bursts would start after the improvements are 
made to NIOBE and the already available data would be of 
very limited use. 

An additional possibility is that similar improvements 
be made to other existing bar detectors. The data corre- 
sponding to GRB times from the five detectors could then 
be co-added and the time for a possible detection could re- 
duce to only 25/\/5 ~ 10 years assuming that all detectors 
cover all GW/GRB events and that they all achieve similar 
characteristics to the improved NIOBE detector. 

Finn, Mohanty & Romano (1999) have independently 
suggested a similar idea to that presented in this article, re- 
garding the future capabilities of the LIGO detectors. They 
suggest that correlating the output from two LIGO detec- 
tors previous to GRB times could result in a 95% confidence 
detection of an association between GRBs and GWs after 
the observation of ~ 1000 GRBs. If no signal is present 
then increased numbers of observed bursts could place lim- 
its on the possible mechanisms for GRBs. Conceivably, a 
detection could be possible within 4-5 years of the begin- 
ning of full LIGO operation. However, if the GW waveform 

and GRBs. The lowest limit on the time lag is 1 s for NIOBE due 
to the ZOP processing algorithm (the sampling time is Is). We 
allow for a further 1 s uncertainty introduced due to the difficulty 
in defining the 'GRB time' as discussed in Sections 2 and 4. 



is known (i.e. the GW/GRB mechanism is determined to 
some degree), a detection of GWs may occur in less time 
since matched filtering could be used to analyse the signal 
more effectively. 

The importance of the above results is clear. If improve- 
ments are made to existing bar detectors, a detection of 
GW/GRB association may be possible in a similar amount 
of time as for twin LIGO-like detectors using the technique 
outlined in Sections 2 and 3. Such detection, whether by bar 
detectors or by laser interferometers, may be the only ob- 
servational technique available to place constraints on GRB 
mechanisms aside from measurements of fading optical GRB 
counterparts. This distinct advantage comes with the loss 
of polarisation information and information about the time 
delay between GW and gamma ray emission. These are, un- 
fortunately, important parameters for distinguishing gamma 
ray burster models (Kochanek & Piran 1993). Also, despite 
the fact that the results in Figs. 3-6 are quite distinct, our 
ability to distinguish the different model parameters is lost 
in the co-addition procedure. However, the importance of 
these losses is not great when we remind ourselves that GWs 
have still not been directly detected; both the lower cost 
improvements to existing bar detectors and the high cost 
development of LIGO detectors are of great importance to 
the study of GRBs. 
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